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ABSTRACT 

On the basis of Borel resummation, we propose a systematical improvement of 
bounce calculus of quantum bubble nucleation rate. We study a metastable super- 
renormalizable field theory, D dimensional O(N) symmetric 4 model (D < 4) 
with an attractive interaction. The validity of our proposal is tested in D — 1 
(quantum mechanics) by using the perturbation series of ground state energy to 
high orders. We also present a result in D = 2, based on an explicit calculation of 
vacuum bubble diagrams to five loop orders. 
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1. Introduction 



In this paper, we propose a new approach for the tunneling phenomenon in 
a metastable super-renormalizable field theory, with aiming at a systematical im- 
provement of the bounce (or instanton) calculus [1,2]. We shall study a D — 1 di- 
mensional system (D < 4) whose Hamiltonian is defined by 



where is an N component real scalar field (0 2 = 0-0) and ir is its conjugate 
momentum. The Hamiltonian (1.1) may be regarded as a Ginsburg-Landau like 
effective theory in which the classical time evolution of order parameter is de- 
termined by the first term. Our present approach therefore might become relevant 
for tunneling problems in condensed matter physics. 

We study an imaginary part of the vacuum energy density. Since the potential 
energy in (1.1) is not bounded from below, the quantum tunneling makes the 
naive ground state (vacuum) = metastable. Therefore, physical quantities 
have to be defined by an analytic continuation from a negative g. In particular, 
the continuation produces an imaginary part of the vacuum energy density Im£, 
which can be regarded as a total decay width per unit volume of the quasi-vacuum, 



The standard method for calculating such a tunneling amplitude in field theory 
is the bounce calculus [1,2]. The existence of the bounce solution, an extreme of 
the Euclidean action with one negative Gaussian eigenvalue, signifies a decay of 
the vacuum due to quantum tunneling [1,2]. The leading bounce approximation 
furthermore gives a quantitatively reliable estimation of the tunneling amplitude 
in the weak coupling region g <C 1. 

| One can find a detailed account on this point in [3]. 





i.e., r = -2lmS/h: 
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The systematical evaluation of higher order corrections to the leading bounce 
result, however, is difficult. One has to consider "interactions" among bounces and 
an integration over the quasi-collective coordinates. Simultaneously, perturbative 
corrections around bounces have to be taken into account. One also has to resolve 
the "mixing" of those two effects. Even if this difficult task could be done, such 
an "instanton expansion" is likely to be an asymptotic expansion, thus we do not 
expect the quantitative validity in the strong coupling region. 

We tackle this issue, i.e., a systematical quantitative improvement of bounce 
calculus, from a completely different viewpoint. Namely we utilize an information 
of the conventional perturbation series around the naive vacuum to evaluate the 
tunneling rate. This should sound strange because usually the quantum tunneling 
is regarded as a non-perturbative phenomenon. However, extensive studies on 
the large order behavior of the perturbation series [4] have revealed an intrinsic 
connection between the quantum tunneling and the nature of perturbation series. 
This connection is the backbone of our approach. Technically, we utilize the Borel 
resummation method [4] for extracting the tunneling rate from the perturbation 
series: Singularities of the Borel transform are expected to reproduce an imaginary 
part of the vacuum energy density. 

In this approach based on the conventional perturbation expansion, every steps 
of calculation (such as the renormalization) are well-understood and, in principle, 
the order of approximation can be systematically increased. To our knowledge, 
this kind of approach to tunneling phenomena in quantum mechanics was initiated 
in [5] on the basis of another kind of resummation method. See [6] and references 
therein. Our present proposal can be regarded as the natural generalization of [6] 
to quantum field theory. 

We will take the "natural unit" k — c — 1 in what follows. 
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2. Bounce Calculus 



Let us first recapitulate the result of bounce calculus [1,2] of the tunneling am- 
plitude. It gives an important information on the nature of a Borel singularity and, 
simultaneously, provides the "standard" with which our results can be compared. 

The Lorentzian action corresponding to Hamiltonian (1.1) is given by 



S[<f>] 



d D x 



\d^<t> - \m 2 <t> 2 + ^(0 2 ) 2 



(2.1) 



where we have defined the Euclidean action Se[4>\ (the time coordinate of Eu- 
clidean spacetime is defined by x = it). We also assume an appropriate counter 
term S coun t.[(f)\ to remove the ultraviolet (UV) divergences: 



S, 



count. 



d D x 



N 
~2 



d D k 



«(2tt) 



D 



ln(m 2 — k 2 — ie) + 



(2.2) 



The structure of abbreviated terms depends on the spacetime dimension D. See 
Eq. (5.1) as the example for D = 2. 

We repeat the procedure in [2] to yield the imaginary part of the vacuum energy 
density for g 1, 



ImS(g) 



Q \ {D+N-l)/2 

= ~A N C D , N ( A) e-Sofi 

bounce \ ZTTg J 



(2.3) 



where dimensionless combinations 8 = S/m and g = g/nrft have been intro- 
duced. 

We briefly explain how the various factors in (2.3) emerge. One expands the Eu- 
clidean functional integral around the spherically symmetric bounce c (r), which is 
a solution of the Euclidean equation of motion, — A0 c (r)+m 2 c (r)— g0 c (r) 3 /3! = 0. 
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The spherically symmetric bounce has the least action [7] whose value is numeri- 
cally given by (S E [(f>c] = S /g)* 



I for D = 1, 

s o = = < 35.10269 for D = 2, 
113.38351 for£> = 3. 



(2.4) 



Gaussian integrations around the bounce except the zero modes give rise to the 
coefficient Cd,n 

C DjN = e ^count.[0c] m -(D+iV-l) (4 _ D) {N-l)/2 

x |det'(-A + m 2 - 9^/2) f 1/2 [det'(-A + m 2 - £0 2 /3!)] ~ (iV - 1)/2 . 

(2.5) 

The determinant can be evaluated analytically in D = 1 and numerically in D = 2 
and Z) = 3 [81: 



C D 



N 



D 



-1/2 



P«(l/3)/4] 



-(iV-l)/2 



2^ x (2 v / 3) jV_1 for D = 1, 

0.3503 x (1.652) iV - 1 for D — 2, 
10.189 x (l^Q)^- 1 for D = 3. 



(2.6) 



Since the quadratic operator — A + m 2 — g4> 2 c /2 has one negative eigenvalue [2,7,8], 
the square root of the eigenvalue produces a factor ±i. A physical requirement 
that we are concerned with a decaying process specified the branch in (2.3). The 



Jacobian from a zero mode to a collective coordinate is a/ So/(2ivg). Because there 
are D + N — 1 collective coordinates (the spacetime position and the direction 
of bounce in N dimensional internal space), the power in (2.3) is resulted. The 
integration over the bounce direction in the internal space gives a factor An, the 
half area of N — 1 dimensional unit sphere, An = tt n ^ 2 /T(N/2). Finally the 



* The numbers, I4, Iq, Dr(1) and Dr(1/3) are defined in [8]: Eq. (58) of [8] should be read 
as At 1 - = (4 - D)D(l/3)/4 for a general D. 
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integration over the position of bounce produces a factor, the system volume times 
the time period, which is divided to give the energy density: This completes our 
quick review of the bounce result (2.3). 

In our approach, we do not need the explicit value of the determinant fac- 
tor (2.6). It will be used merely for a comparison. On the other hand, the value 
of the bounce action (2.4) and the number of collective coordinates are important 
because they tell the position and the strength of the nearest Borel singularity to 
the origin of Borel plane. They can be easily obtained: One only has to solve a 
one dimensional differential equation to find the bounce action, and a symmetrical 
consideration fixes the number of corrective coordinates. 



3. Resummation of Vacuum Bubbles 

The leading bounce result (2.3) is reliable for the weak coupling region g <^1. 
In this section, we present our proposal which is expected to work even in the 
strong coupling region. 

We start with the conventional perturbative expansion of the vacuum energy 
density, namely a sum of the vacuum bubble diagrams, 

oo 

£(«?)- J>S?\ (3.1) 

n=0 

where we have assumed an appropriate renormalization which makes c n s finite. 
For super-renormalizable cases D < 4, only first several CfiS £1X6 UV divergent. 

From the perturbation series (3.1), we construct the Borel (more precisely 
Borel-Leroy) transform: 

oo 

The argument of gamma function in the denominator has been chosen so that 

★ However, if one has perturbative coefficients to sufficiently higher orders, even those infor- 
mation from bounce calculus may be deduced from the perturbation scries. Sec [6] 
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a singularity of the Borel transform nearest to the origin becomes a square root 
branch point (see (3.5)). 

When the coupling constant g is negative in (1.1) and (2.1), namely with a 
bounded potential, the proof of the Borel summability of Green's functions [9] 
may be generalized to the vacuum energy density. With our definition (3.2), the 
Borel summability is expressed as 

oo 

^ = h)» / <^ (D+ " )/2 - lB (-4 (3-3) 

Now, we assume that the Borel transform B(z) is analytic on the complex z plane 
and possible singularities of B(z) exist only on the positive real axis. Precisely 
the same assumption on the Borel transform of Green's functions has been made 
in the perturbation approach to the critical phenomena [11,4]. We also assume 
the behavior of B(z) at z = oo is moderate enough and the convergence of the 
integral (3.3) is determined by the exponential factor. 

Under these assumptions, we may rotate the integration contour in (3.3) as 
z — > e^ m z as it goes along the upper (lower) side of the positive real axis. This 
operation gives rise to the analytic continuation for g > 0, 

oo 

£(<?) = ^v)72 / dze-^z^ D+ ^ 2 - 1 B(z±ie), (3.4) 
o 

which is our basic equation. 

With (3.4), the leading bounce result (2.3) of the imaginary part for j C 1 
implies that there exists a singularity of Borel transform at z = So [4]. Because of 
the structure of (3.4), it is the nearest singularity to the origin. A comparison of 
(2.3) with (3.4) shows that the singularity is a square root branch point: 

1 9 1/2 
B(z) = -—A N C DiN {2 ^J +N _ 1)/2 (S - z)-V* + ■■■. (3.5) 

When substituted in (3.4), the branch cut (3.5) reproduces the imaginary part (2.3) 
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for g <^ 1. The sign of the singularity (3.5) itself is not determined solely from 
the comparison of (2.3) with (3.4) because there are two possible choices of the 
integration contour (upper side or lower side of the positive real axis). As a consis- 
tency test, we note that the expansion of (3.5) with respect to z yields negative c n s 
in view of (3.2), which is nothing but the large order behavior of the perturbation 
series [4]. As we will see below, this is a correct property of actual perturbation 
coefficients when we follow the — ie prescription. This consideration also shows 
that we have to take the upper contour to obtain a negative imaginary part, which 
is the physical one. 

Basically we construct the Borel transform (3.2) from the actual perturbation 
series c n and substitute it in (3.4) for extracting the "non-perturbative" informa- 
tion. However the radius of convergence of the series (3.2) is finite (= So) due to 
the singularity (3.5) and we have to analytically continue the series (3.2) outside 
the convergence circle to perform the Borel integration (3.4). This is an impossible 
task without knowing all the perturbative coefficients. However, as is well-known, 
this difficulty can be avoided by the conformal mapping technique [10]. We thus 
introduce a new variable A by 

i = 4S »(TTaF (3 ' 6) 

The point is that the convergence circle of the series (3.2) in terms of A is now 
a unit circle, within which the whole cut z plane is mapped. In particular, the 
real axis z > So is mapped on the circle |A| = 1. Therefore we may use a finite 
order truncation of the series of A in the Borel integration (3.4) order by order. 
Then we may expect the sequence converges to (3.4) (as was the case in quantum 
mechanics [6]). 

In terms of A, the Borel transform (3.2) is expressed as 
Combining (3.4), (3.6), and (3.7), and parametrizing the unit circle by A = e l °, we 
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find the Pth order approximation of the imaginary part, 



'SoV D+N)/2 j M ( So 1 \ sinfl/2 f . (3-8) 
T ) r 9 6XP oo^ ) cos^+i 6>/2 g 4 Sm ^ 



Note that this is solely expressed by the first P perturbative coefficients c n and 
the value of bounce action So- 

Our formula (3.8) relates the value of vacuum bubble diagrams c n and an 
imaginary part of the vacuum energy density In quantum field theory, usually one 
never seriously compute the vacuum bubble diagrams because one should be free 
to adjust the origin of a real part of the energy density and thus it has no direct 
physical meaning. On the other hand, the imaginary part of the vacuum energy 
density has a definite physical meaning as the decay width of the quasi-vacuum. 
Therefore the latter should not be modified under a change of convention for the 
former. How can this general consideration and the formula (3.8), which relates 
those two quantities, be reconciled? 

In fact, there is no real contradiction. For example, suppose that we perform 
a finite renormalization of the vacuum energy density to eliminate the first Q per- 
turbation coefficients. (Note that we cannot take Q = oo because a simple sum 
of the perturbation series is diverging.) Then the new Borel transform would be 
defined by 

oo Q 

Note that the last finite series (i.e., polynomial) does not develop a singularity and 
thus cannot contribute to any imaginary part via the Borel integral. Therefore the 
imaginary part, as it should be, is invariant under such a change of the origin of 
the energy density. Contrary to its peculiar looking, the formula (3.8) is workable 
in this way. 
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4. Quantum Mechanics 



We can extensively test the validity of our master formula (3.8) in quantum me- 
chanics D — 1. In this case, perturbative coefficients of the vacuum energy, namely 
the ground state energy, to very high orders are available. The most efficient way 
for computing them is the recursion formula method [12]. By generalizing it to 
the O(N) symmetric model, we have computed c n to n = 50. The first several 
coefficients are* 



N N(N + 2) N(N + 2)(2iV + 5) 

°° ~ "2 ' Cl ~ 96 ' ° 2 ~ 4608 
N(N + 2)(8N 2 + A3N + 60) 

C3 " 



221184 f4 . 

iV(iV + 2) (168iV 3 + 1437iV 2 + 4270iV + 4420) v ' ; 

C4 = 



C5 = 



42467328 

N(N + 2)(1024iV 4 + 12277iV 3 + 57668iV 2 + 126128iV + 108480) 

2038431744 ' 



The exact complex quasi-ground state energy is also available by a numerical 
diagonalization of the Hamiltonian in a (unconventional) Hilbert space with a 
rotated boundary condition (see, for example [6]). Therefore we can compare our 
formula (3.8) and the bounce result (2.3) with the exact value of the imaginary 
part. 

In Fig. 1, we have plotted the result of (3.8) with N = 1 and P = 4, P = 5, 
and P = 15, respectively.^ The solid line is the exact value obtained by the numer- 
ical diagonalization of the Hamiltonian. The imaginary part is normalized by the 
leading bounce result (2.3). The broken line, which is depicted for a comparison, 



* It is worthwhile to note that all the coefficients (except Co) are proportional to N + 2; 
this holds in an arbitrary dimension D. This can be proven by noting that the partition 
function for N = — 2 can be expressed by a single component fermionic system with a trivial 
interaction (ipip) 2 = 0. 

f This result has already been reported in [6] . 
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is the bounce result including two loop radiative corrections [13], 
lm£(g) 

J bounce plus two loops 

^L) N, \- S ,H ( 1 _ ^ + 5W + 20- N <4 ' 2) 



= - AnC ^\M) ^'•' S { 1 576 », 

We see an excellent convergence of (3.8) to the exact value and, as was announced, 
the proposal in fact gives rise to the improvement of bounce calculus. The plotted 
range of the coupling constant was determined by a criterion that the real part 
of ground state energy (which can also be compute numerically) is lower than the 
potential barrier. Therefore the plotted range may be regarded as the quantum 
tunneling (not classically sliding) region. The agreement with the exact value is 
better for g larger but this is physically reasonable: When g becomes smaller, the 
potential barrier becomes higher and wider and thus the tunneling phenomena is 
difficult to be detected from a perturbative expansion around the potential origin. 

Fig. 2 is the same as Fig. 1, but for N = 2. In this case, P = 15 shows a small 
excess in the weak coupling region. However we verified that (3.8) eventually 
converges to the exact value, as indicated by the plots of P = 30 (squares). 

In both figures, the bounce calculus including two-loop corrections (the broken 
line) is giving a rather nice fitting of the exact value. Therefore one may wonder 
whether the imaginary part is almost saturated by the perturbation corrections 
around the bounce, or, the multi-bounce contribution is crucial in this coupling 
constant region. It is thus of interest to study how large the multi-bounce contribu- 
tion is. We may estimate the two bounce contribution by studying the interaction 
between bounces and the quasi-collective coordinate (e.g., separation between two 
bounces) integration [14]. The interaction is attractive when they have the same 
orientation and thus the two bounce contribution dominates the functional inte- 
gral when the separation is small. However, since the notion of multi-bounce is 
meaningful only for a large separation, one may define the partition function by 
an analytic continuation from the negative g [14]. In this way, we find for N — 1 



| The sign of the contribution cannot be fixed by this prescription alone. 
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lm£(g) 



. two bounces 



±2 [ Im ^) 



. bounce 




(4.3) 
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which is 42% of the one bounce contribution for g — 4, the boundary of the 
tunneling region. Therefore, if one pursues the bounce calculus, the inclusion of 
multi-bounce contributions should be crucial in this coupling constant region. 

From the excellent agreements in Figs. 1 and 2, therefore we conclude that 
our approach based on the perturbation series around the trivial vacuum is in 
fact taking into account the multi-bounce contributions. In other words, our as- 
sumption (3.4) contains a statement that an imaginary part of the vacuum energy 
density is completely saturated by Borel singularities of the perturbation series 
around the trivial configuration. 

As the final application of the quantum mechanical model, we may test our 
assertion made in (3.9) that a finite subtraction of the perturbation series does not 
affect the imaginary part when P is sufficiently large. Fig. 3 is the same as Fig. 1 
but the first two coefficients cq and c\ are set to zero by hand. We see that the same 
kind of convergence behavior even with this subtraction (except the intermediate 
oscillating behavior; note that P = 4 is better than P = 5) and P = 15 reproduces 
almost the same result as Fig. 1. 

5. Tunneling on Line 

After observing our proposal works quite well in quantum mechanics, let us 
try to apply our formula (3.8) to D — 2, that represents a certain one dimensional 
system (line or wire). For D = 2, we assume the following counter terms in (2.2): 




x {f/|^ b ( m2 -* 2 - ! ')-> + 2 »»/ 

+ ai NiN + 2) «\J <(ar)»m»-*»-fe }' 



i(2n) 2 m 2 — k 2 — it 



d 2 k 1 



(5.1) 
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where the first and third terms remove the zero-point energy to two loop orders, and 

c 

the second is the counter term for the one loop self energy. Since the present system 
is super- renormalizable, the counter terms (5.1) remove all the UV divergences to 
higher orders. 

To find the lower order perturbative coefficients c n in (3.1), we have explicitly 
calculated the vacuum bubble diagrams to five loop orders. With the counter 
terms (5.1), cq = c\ = 0, and there are five diagrams to be evaluated (Fig. 4). 

We have the following numbers: The diagram (a), 



c 2 = K — x 8.833 895 1 x 1(T 5 . (5.2) 

o 



_ = _ iV ( iV + 2 ) x 8.833 8951 x 1(T 5 . 

The diagram (b), 



N(N + 2)(N + 8) a , , 

c 3 = 1 — --21 — L_L x 3.012 767294 x 1(T 6 . (5.3) 



The diagrams (c), (d), and (e), respectively, 



N(N + 2)(N 2 + 6iV + 20) 

c 4 = ^— x 5.657478 653 058 x 10" 

81 



N{N + 2) 2 



x 1.006 825 50 x 10" 7 (5-4) 



_ V + 2)(57V + 22) 
81 

In the above numbers, we have verified only the last digit contains the roundoff 
error due to numerical integrations. In the diagram (e), we used the analytical 
expression of one loop triangle diagram evaluated by [15]. The numerical error 
for the diagram (e) is much larger than others: However we verified that the final 



§ Therefore, our coupling constant g is not the physical four point scattering amplitude. The 
finite renormalization, which is required to translate the result into the one in terms of 
the physical coupling constant, might be performed by utilizing the resummed f3 function 
in [11]. 
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result in Fig. 5 is quite insensitive to the precise value of (e) by varying it within 
the estimated error. As the general property, all the diagrams have an overall 
minus sign according to the — ie prescription. Notice that we are computing the 
energy density, instead of the effective action. 

We could not prepare the exact value of the energy density, although this task 
might be done by invoking a lattice formulation in a complex rotated functional 
space. Therefore it is not clear which range of g is the tunneling region in this 
case. As an order of magnitude estimation, we may consider the effective potential 
to two loops which behaves around the origin as 



i-i(iV + 2)x6.18xlO- 4 ^ 2 



m 2 2 + • • • , (5.5) 



where the counter terms (5.1) have been taken into account. To the two loop 
order, the real part of the vacuum energy density vanishes, cq = c\ = 0. Therefore 
we may regard the value of g which changes the sign of curvature of the effective 
potential at the origin as a measure of the tunneling region. This rough estimation 



indicates < g & a/3/ (N + 2) x 30 is the tunneling region. 

In Fig. 5, we have plotted an imaginary part of the vacuum energy density (3.8) 
with P = 2, P = 3, and P = 4. Again it is normalized by the bounce result (2.3). 
Unfortunately, it seems impossible to draw a definite conclusion from Fig. 5: We 
do not observe the convergence behavior and in some region even an unphysical 
result, the positive imaginary part, can be found. It is not clear whether this is 
due to the luck of orders of the perturbation series, or there exists a fundamental 
obstruction for our approach we did not encounter in quantum mechanics. (The 
nature of the nearest Borel singularity in the present case, D = 2 and N = 1, is the 
same as that of N = 2 quantum mechanics, Fig. 2.) If we nevertheless take Fig. 5 
at its face value, it is suggesting that the true tunneling amplitude is much larger 
than the leading bounce result in the strong coupling region. However, this would 
remain a speculation without having much higher order perturbation coefficients, 
which will clarify the real convergence property of our proposal. 
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6. Conclusion 



In this paper, we have proposed a new approach to the tunneling phenomena 
"the decay of the false vacuum" in a super-renormalizable field theory. Our ap- 
proach utilizes the information of conventional perturbation series around the naive 
vacuum = 0. We have verified numerically that the singularities of the Borel 
transform reproduce an accurate tunneling rate in quantum mechanics (D = 1). 
We have also presented a result in D = 2, based on an explicit calculation of vac- 
uum bubble diagrams to five loop orders. Unfortunately the number of orders of 
the perturbation series is not sufficient to make definite conclusions on the con- 
vergence property and the true tunneling rate in this case. Only the higher order 
calculation will answer these questions. 

Although we have only presented results for D = 1 and D = 2 in this paper, 
the formula (3.8) is waiting for the straightforward application to D = 3 = 1 + 
2 dimensional system (tunneling on plane). We expect a calculation of perturbative 
coefficients to six loops is tractable because the analytical structure of one loop 
diagrams in D = 3 is rather simpler than that of D = 2 [15]. As a consequence, 
definite conclusions on a convergence property of the method and on the true 
tunneling rate may be drawn. We hope to come back this problem in a near 
future. 

Another possible test of our proposal is offered by an interesting field theoretical 
model in [17], for which both the bounce calculus and the perturbative calculation 
to very high orders are possible [18]. 

Finally we briefly comment on a possible generalization to the just renormal- 
izable case, D = 4 = 1 + 3. The first trouble one encounters is that all the 
vacuum bubble diagrams have an overall (i.e., not sub-diagram) UV divergence. 
A simple renormalization which sets all the coefficients c n zero order by order, 
appears meaningless* Therefore a natural procedure is that one first introduces an 

* A similar situation occurs when one evaluates an imaginary part of the vacuum energy 
density due to the UV renormalon. See the third reference of [16]. 
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UV cutoff while hoping the imaginary part (3.8) is finite and independent of the 
regularization for P — > oo. Next we are worried about an emergence of the UV 
renormalon [16], another known source of the Borel singularity. Interestingly, it 
does not emerge in our present case because 4 model with a negative g is asymp- 
totically free (no Landau pole)? Of course, since 4 model in D = 4 is trivial when 
the UV cutoff is removed, the relation to the "true" imaginary part is not obvi- 
ous. Presumably, our proposal for the four dimensional 4 model is meaningful 
only with an UV cutoff, but such an UV cutoff is naturally provided in condensed 
matter physics. 

The work of H.S. is supported in part by the Ministry of Education Grant-in- 
Aid for Scientific Research, Nos. 08240207, 08640347, and 08640348. 
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FIGURE CAPTIONS 

1) Ratio of an imaginary part of the quasi-ground state energy to the leading 
bounce result in D = 1 and N = 1 case. The solid line is the exact numerical 
value. Eq. (3.8) with P = 4 (circle), P = 5 (filled square), and P = 15 (filled 
circle), is plotted. The broken line is the leading bounce result with the two 
loop correction (4.2). 

2) Same as Fig. 1, but for N = 2. Eq. (3.8) with P = 30 (square) is also plotted 
for a comparison. 

3) Same as Fig. 1, but with the first two perturbation coefficients are set by 
hand to zero, co = c\ = 0. 

4) Vacuum bubble diagrams with the counter terms (5.1) to five loop orders. 

5) Ratio of the imaginary part of the vacuum energy density computed by (3.8) 
to the leading bounce result in D = 2 and N = 1 case. P = 2 (circle), P = 3 
(filled square), and P = 4 (filled circle), are plotted. 
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